Landau Levels and Quasiparticle Spectrum of Extreme Type-II Superconductors 
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At high magnetic fields and low temperatures, numerous 
extreme type-II superconductors exhibit Landau quantization 
of electronic motion. We present an analytic construction of 
the quasiparticle spectrum in this regime, based on the high- 
field expansion. The spectrum is gapless and is separated from 
the familiar low-field regime by a quantum level-crossing tran- 
sition. Such low energy excitations should lead to observable 
effects in thermodynamics and transport. 

Recent observations of the de Haas-van Alphen 
(dHvA) oscillations in many extreme type-II supercon- 
ducting materials! clearly indicate the presence of Lan- 
dau level (LL) quantization within a high-magnetic field 
(H), low-temperature (T) "pocket", surrounding the 
H C 2(0) point in the H — T phase diagram (Fig. 1). The 
size of this pocket in conventional type-II systems (like 
Nb) is expected to be negligible. This is so because, 
within the BCS theory, the scale of the cyclotron splitting 
between LLs near H c2 (0), w c ~ w c2 (0) = eH c2 (0)/m*c, is 
set by the condensation energy, ~ T 2 /Ep, and should be 
much smaller than either the thermal smearing, ~ T, or 
the gap A(T, H), the scale for both being T c q. Additional 
smearing due to disorder, T, makes this tiny pocket effec- 
tively irrelevant. In clean, intrinsically extreme type-II 
systems, the situation is significantly different. We de- 
fine the extreme type-II systems as those with the slops 
of the upper critical field at T c0 , H' c2 > 0.3 Tesla/Kelvin.El 
In such systems w C 2(0) is comparable to T c q and there is 
large pocket in the phase diagram in which the LL struc- 
ture within the superconducting phase is well defined, 
i.e. uj c > A(T, H),T, T. Numerous superconductors be- 
long to this extreme type-II family: high-T c cuprates, 
A15's, boro-carbides, many organics, etc. The bound- 
aries of the "extreme" pocket, H* and T*, defined by 
u>c ~ A(T = 0, H) and u c ~ T, extend to H as low as 
H*~ 50%H c2 (0) and T as high as T* ~ 30%T c0 (Fig. 
l).El 

Within this "extreme" pocket the solution to the su- 
perconducting problem must fully incorporate the LL 
structure of the normal state. This leads to a set 
of Bogoliubov-de Gennes (BdG) equations for the LL 
quantized quasiparticles in presence of the gap-function 
A(r), describing the Abrikosov vortex latticema There 
are two main sources of difficulty in solving these equa- 
tions. First, the number of LLs involved in the pair- 
ing is typically 40-100 and the matrix elements of the 
gap function between these states, A nm (q), are in gen- 
eral quite complex. A more serious difficulty, however, 
is that the basis which diagonalizes the BCS Hamilto- 
nian involves combined "rotations" both in the Nambu 



space and the LL basis. To illustrate why this fact seri- 
ously impedes analytic progress we first note that, to the 
leading order in A/w c <C 1, the quasiparticle excitation 
spectrum near the Fermi surface (FSA can be found ex- 
actly: E = ±y / e n (k z ) 2 + |A, m (q)| 2 .B It is tempting to 
generalize this and conclude that the full solution takes 
the form: E = ± yEf^ + \D nii \ 2 , where T, nfi (q, k z ) and 
Dnn(c[,kz) are the normal and pairing self-energies, re- 
spectively. The physical meaning of such solution would 
be that we first rotate the LL basis, {n} — > {fi}, and 
then diagonally pair up these new "normal" quasiparti- 
cles. This would lead to a simple description, suitable 
for applications of our standard theoretical machinery. 
Unfortunately, this conclusion would be incorrect: be- 
yond the leading order, the full solution cannot be put 
in the above simple form as the "normal" and "pair^ 
ing" self-energies cannot be simultaneously diagonalized.tj 
This leads to a loss of physical transparency and makes 
it difficult to interpret results in familiar terms. Mc 
the work beyond leading order has been numerical.l 

In this Letter we introduce an analytic approach which 
allows for a transparent and systematic evaluation of cor- 
rections to the leading order results. The basic idea can 
be illustrated by an example of an interacting Fermi liq- 
uid: there various properties are in general very diffi- 
cult to evaluate. However, as long as we are interested 
in the low T thermodynamics and low energy transport, 
these properties can be accurately described by the quasi- 
particle representation of the excitations near the FS. 
This leads to the low T thermodynamics which is in 
one-to-one correspondence with that of non-interacting 
fermions, but with a modified dispersion relation and a 
reduced spectral weight. In this paper, we construct such 
a "quasiparticle" representation of the present problem. 
Of course, there are no interactions in our case but our 
aim is similar. The original partition function is replaced 
by the one of "quasiparticles" whose energies do have 

the desired form E = ± \J s 2 ln + Z 2 1 A „„ | 2 , and which 
faithfully represents the low T (or low energy) properties 
of the original problem. Here e nn (q, k z ) and A nn (q, k z ) 
are the renormalized normal and pairing self-energies and 
Z is the "quasiparticle" renormalization factor. All can 
be computed to the desired accuracy in A/uj c from a 
general expression which we provide. Using this "quasi- 
particle" representation we derive the following results: 
first, we evaluate the leading order corrections to "nor- 
mal" and "pairing" self-energies and find that they are 
of order ~ A 2 /oj c and ~ A 3 /w 2 , respectively. Therefore, 
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for A/uj c <C 1, the features of the spectrum obtained in 
the leading order are only moderately affected. In partic- 
ular, we show that the quasiparticle excitations remain 
gapless to all orders in perturbation theory until A/ui c 
exceeds certain critical value. As A/uj c increases further 
(H < H*), the original LL structure is obliterated by 
strong LL mixing induced by off-diagonal A nm 's. This 
exact result shows that the gapless "extreme" pocket is 
truly a novel state of type-II superconductors and is sep- 
arated from the gapped low-field mixed phase by a quan- 
tum level crossing transition. In this context, we point 
to an important difference between 3D and 2D systems. 
We next present our numerical results for the excitation 
spectrum and find that they are in agreement with the 
analytic results. We hope to inspire experiments which 
could test our results and provide novel insight into the 
physics of high-field, extreme type-II superconductors. 
Our starting point are the BdG equations: 

m 

E V k z ,q,n = — £ n{k z ) v k z ,q,n + ^intifll""^^ ) (1) 

m 

where e n (k z ) = h 2 k 2 /2m + nhu c — fj,, n is the LL index, 
huj c /2 has been absorbed into the chemical potential, /i, 
q is the magnetic wavevector and k z is the wavevector 
along the field direction^ A„ m (q) =J2j A£, m (q) is the 
matrix element of A(r) = Y^j Aj(rju between electronic 
states (fe 2 ,q, n) and (— k z , — q, m).tl N is the number of 
LLs which participate in pairing. The solution of the 
above equations gives the quasiparticle states and the 
energy spectrum within the superconducting phase. 

As one crosses into the superconducting state and A 
becomes finite, there are both perturbative and non- 
perturbative effects in Eq. (Q). We first isolate the 
non-perturbative part by noticing that the particle bands 
(it-bands) with the LL index ri\ = n + p and the 
hole bands (u-bands) with the LL index ni = n — p 
(p = 0,±1,±2,...) cross at energies phio c whenever 
k z = ikpn = y2m(/j — nhuj c ), where kp n are the 
Fermi wavevectors of different quasi one-dimensional LL 
bands. Additional band crossings t ake place at ener- 
gies (p + l/2)huj c and for k z — ±y / k' Fn + mhuj c . At 
the band crossings near k z — ±fej?„ the matrix elements 
A„+ p . Tl _p(q) act non-perturbatively and must be treated 
by a degenerate perturbation theory. The same is true for 
(q) near k z = ±y / k 2 ?n + mhui c . The effect of 
all other matrix elements, A n+p . n (-i)- P ' (q) {p ^ p'), is 
purely perturbative and can be computed by an expan- 
sion in A/oj c . Near H C 2, where A/oj c <C 1, that effect is 
small and we can obtain the analytic form for the BCS 
excitation spectrum near k z = ikp n : 



E n . p = phw c ± ye n (k z ) 2 + |A„ +P ,„_p(q)| 2 , (2) 

where p — 0, ±1,±2, .... A similar expr ession can be 
written for the spectrum around k z ~ ±yJk 2 Fn + mfiuj c . 



Typically, we are interested in the quasiparticle spectrum 
when we discuss low T (or low energy) properties, i.e. 
(T,uj) <C A(T,H) < oj c . In that case only the quasi- 
particles near the FS are important, k z ~ ±kp n and 
p = 0, and it suffices to consider only the E n , p== Q — 
±y/e n (k z ) 2 + \A nn (q)\ 2 bands. The excitations from 
other, p ^ 0, bands (|2j) are gapped by energies which are 
too high, E ntPi £o(k z ~ ±kFn) ~ phu> c and their contri- 
bution to the quasiparticle thermodynamics is negligible 
in the low T (T <C A(T, H) < w c ) regime. This is the 
content of the so-called "diagonal approximation" (DA) 
for the excitation spectrum.tffl 

The main issue addressed in this paper is how to go 
beyond Eq. (|^) analytically. As one moves further below 
H C 2, A/ujc increases (while still <C 1) and the correc- 
tions coming from A„ +Pi „_ p '(q) (p ^ p') must be in- 
cluded. A pedestrian perturbation theory, starting from 
(||), leads to expressions whose physical content is diffi- 
cult to analyze. The reason is that (^) already involves a 
full Nambu rotation. Consequently, all other terms gen- 
erated by the pedestrian perturbation theory are fully 
Nambu rotated. This makes it impossible to disentangle 
the "normal" from the "paring" part of the problem.tl 
In contrast, our standard theoretical machinery works 
best when we can neatly separate the normal from the 
pairing part. We propose here a different method of com- 
puting the excitation spectrum. Consider the partition 
function specified by the solution of BdG equations (ffl) , 
Z = rL, q ,fc, z ( iuJ > q> k z ), where 

Z = J] / dt/j n d'4> n exp{y^ ^njiuj'I - H) nm ^ m } . (3) 

Here the set of Grassman variables {u n , v„}(«w, q, k z ) 
representing the particle (hole) bands of O) has been ar- 
ranged in spinors %p n = [u n v n ). I and 7i(q, k z ) are two 
2N x 2N matrices which are, respectively, a unit matrix 
and a Hamiltonian whose form is evident from the right 
hand side of Eq. ([!]). {lu} is the set of fermionic Matsub- 
ara frequencies. Note that Z(iu>, q, k z ) is a determinant 
of a matrix, Z = dct ||icjZ— H\\. When analytically con- 
tinued, E — > iu>, this is precisely the secular determinant 
of the BdG equations (|l|). Thus, the solutions for quasi- 
particle energies in (Q) are the zeroes oiZ{iu> — > E, q, k z ). 

Within our "extreme" pocket (T < A(T,H) < u c ) 
the main contributions to Z come from k z ~ ikp n and 
E n ,p=o bands. In order to evaluate Z(itu,q, k z ~ kp n ) 
we first integrate out all tp njP ^o(iLL!,q, k z ~ kp n ) in (|^): 

Z = Z pz£0 /#n,o#n,oexp{'0 n>o (iwX-'H)ooV'n,o}- (4) 

Z p ^o{iui, q, k z ~ kp n ) is what one gets from (^) with 
the E n ^ p= Q band excluded. It describes the contribution 
to thermodynamics arising from all other bands, -E„ iP ^o> 
by themselves. Since p ^ bands are far from the FS, 
-E'n,p#o(q, k z ~ kp n ) ~ pTiuJc, such contributions are ac- 
tivated with a large gap, i.e. ~ exp(— pfrui c /T) <c 1, and 
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their contribution is very small. Thus, the dominant con- 
tribution comes from det || (iul— H)qo || ■ Here T o is a 2 x 2 
unit matrix while Hoo = H°+Tl' is a 2 x 2 Nambu matrix. 
The diagonal elements of {Tt°, H'}(iu), q, k z ~ kp n ) are 
{±e n (k z ), ±T, nn (iuj, q, k z )} while the off-diagonal ones 
are {A„„(q), Z?„„(iw, q, k z )} and their complex conju- 
gates (c.c). By itself, TC° determines the spectrum of 
quasiparticles at the FS (described by ip n ,p=o) in the DA 
H! is a correction arising from coupling to other 
spinor bands, i/) n ,p^0i with S nn and D nn being the nor- 
mal and pairing self-energies, respectively. Since the 
p / bands are separated from the FS by phtu c , these 
self-energies contain denominators ~ phuj c and are small 
for A/uj c <C 1. The explicit form is obtained from: 



p> 



(5) 



Ti.0p is a 2 x 2 Nambu matrix whose diagonal elements 
are zero and off-diagonal ones A„ in+ p(q) and its c.c. The 
resolvent appearing in (||) is that of the BdG system (|l|) 
for k z ~ kp n but with the p = spinor band (||) excluded. 

Up to this point there are no approximations and the 
expression for TC is formally exact. The problem is that 
the resolvent for the p ^ bands involves inversion of a 
large matrix. However, for A <C w c , the resolvent can be 
calculated perturbatively, to a desired accuracy in A/lu c . 
In general, it is useful to first perform the Nambu rotation 
for the p 7^ bands to get rid of the degeneracies between 
the n + p particle and n — p hole bands and then proceed 
with the perturbative expansion. It turns out that this 
is not necessary if we are only after the leading order 
corrections. In this way we obtain: 



E 

p#0 



P (q) 



iuj + e 



n+p 



&nn 



A n,n+p(q)A; + p in+p ,(-q)A n+ p/ |n (q) ^ 



p^0,p'#0 



{iuj + e n+p (k z ))(iuj - s n+p '(k z )) 



^nn{i<jJ, q, k z ) and D nn (iuj, q, k z ) are the normal and 
pairing self-energies, respectively. They have a frequency 
dependence arising from the coupling to the p ^ spinor 
bands through the matrix elements A„ in+J) (q). We can 
view them as self-energies of some "interacting" system 
and proceed to construct the "quasiparticle" representa- 
tion of (|3j) . This construction uses a Fermi surface expan- 
sion of (0). First, we determine the renormalized "quasi- 
particle" energies from e„„(q, k z ) — e n (k z ) + Y, nn (iuj — 
e„„(q, k z ), q, k z ). To the leading order in A/uj c this gives: 



£rm(l) k z ) PS £n(k z ) + 



|A,; : „ + p(q) 

p#0 



piVc 



(7) 



We also evaluate the "quasiparticle" renormalization fac- 
tor, Z^ 1 = 1 — dT, nn /diuj\i u —g nn ; to the leading order 
z n » 1 - J2 P ^o \&n,n+p(q)\ 2 /p 2 ^c- Similarly, we can 



compute the "renormalized" shape of the FS. The "renor- 
malized" Fermi velocity is given by: 



VFn(q) 



1 



(q)l 



V F „ 



p^O 



P 2 uj1 



(8) 



s nn and VFn determine "quasiparticle" density of states 
(DOS) at the FS and their q-dependence describes the 
"broadening" of LLs by the potential scattering induced 
through the spatial non- uniformity of A(r). This effect 
is secondary to the leading pairing- induced "splitting" of 
LLs in Eq. (||). Finally, the renormalized paring matrix 
element at the FS is: 



E 



A n „(q) = A„„(q)- 
A„,„+p (q) A; +p ^ +p , (-q) A n+P > } „ (q) 



pp'ul 



(9) 



This "quasiparticle" construction is not an ordinary 
perturbation theory. It keeps the "normal" and "paring" 
self-energies separate by postponing the Nambu rotation 
until the last step. It can be viewed as a "quasiparticle" 
representation of the BSC quasiparticle spectrum. Such 
representation is valid as long as we are interested in 
excitations near the FS. 

The above construction can be carried out to any de- 
sired degree of accuracy in A/lo c , although the amount 
of algebra rapidly increases. Certain general statements, 
however, can be made. In particular, the excitation spec- 
trum determined by the BdG equations (|l|) is gapless in a 
finite interval A/lj c 6 [0, x c ], where x c ^ is some criti- 
cal value. This property can be demonstrated by noticing 
first that all A„ m (q) have a set of common zeroes whose 
location in the q-space depends only on whether n + m 
is eves or odd. These are so-called Eilenberger zeroes, 
{qf }.H In the leading order (§) the gap is determined by 
the diagonal matrix elements, A rm (q), and it vanishes at 
q = qf as well as at numerous other points {q™™} within 
the magnetic Brillouin zone (MBZ). As the matrix ele- 
ments coupling different VVp' s are turned on, it is clear 
from (||) that D mi (q = qf ) involves only terms with an 



odd number of matrix elements A 



n+p,n+p' 



(q), i.e. 



Dnn ^ ^n,n-\-p^-n-\-p,n-\-pi ^n+pi ,n+P2 



A 



Tt+pfc,Tt 



(10) 



Note now that it is impossible to form such a string ( JTc| ) 
of odd number of A's without at least one A„ + p/ in+p » (q) 
having n + p' + n + p" even. But such a matrix ele- 
ment also must vanish at same {qf } as A nn (q). Conse- 
quently, at these points, q = qf , the spectrum remains 
gapless. This statement is correct to all orders in pertur- 
bation theory and therefore is exact as long as pertur- 
bative expansion itself is well-defined. We expect that 
the radius of convergence of the perturbative |[ £ ; xpansion 
extends from A/u> c — to A/iv c = x c ~ lJafca At this 
critical point, the perturbation theory breaks down due 
to band crossings between neighboring LL branches. For 
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A/uj c > x c the gaps will open up at {qf } signaling the 
destruction of the LL structure by strong mixing due to 
increasing A. Ultimately, for A/w c 3> x c , there is a 
crossover to the low- if regime of mini-gapped states in 
well-separated vortex coresa, as established by Norman 
et al.u. This phenomenon of band crossings and disap- 
pearance of gapless excitations has been already noticed 
in numerical solutions of the BdG equations.EI Here we 
elucidate its physical origin through analytic arguments. 

We now compare our analytic results with the numer- 
ical solution of the BdG equations (Q). Such solution is 
very time consuming in the 3D case so we have solved (|l|) 
for the 2D case. Our analytic results can be adapted to 
the 2D case but with one key difference: in a 2D system 
the FS is a line not a sheet. Consequently, the gapless 
and "near" gapless points in the MBZ might not be lo- 
cated on the FS. In contrast, the Fermi "sheet" of a LL 
quantized 3D material contains all of the MBZ. As a re- 
sult, the gapless character of the excitation spectrum will 
be less pronounced in 2D systems even though the effects 
of LL quantization itself are more in evidence there than 
in quasi 3D materials. This difference between 3D and 
2D systema.jshould lead to observable effects in dHvA 
oscillationsoaO, thermodynamics and transport. 

We obtain the spectrum as a function of A for differ- 
ent n c = fjb/1uo c . The eigenvalues are grouped in pairs 
centered around the LL energies. As A grows, the bands 
start to broaden until eventually the gap between them 
goes to zero and the band-crossing occurs. The sequential 
eigenvalues have a mirror-like shape with the neighboring 
bands. We fixed \x such that {qf } points remain gapless. 
In 2D, if is adjusted so that an originally gapless point 
in the MBZ (^) is on the FS, then the gap opened by the 
off-diagonal terms should be due to the pairing term and 
thus ~ A . In general, however, the normal contributions 
at q = {q™"| gapless points are finite and their gaps will 
be ~ A 2 . A fit of the normal and pairing self-energy 
contributions for q = q™" as a function of /i, shows that 
there is a value of /i for which the coefficient of the A 2 
term vanishes and therefore close to this point the A 3 
term prevails. This provides a criterion-Jpr locating the 
"renormalized" 3D FS at all {qf ,q^ n }B 

An effective measure of the importance of low energy 
excitations is the density of states (DOS). The shape of 
the DOS at low energies is determined by .the first and 
second order zeroes found in the spectrum J13 The exact 
DOS for small A and DOS within the DA are basically 
equal since A/cj c <C 1. There is therefore a large DOS at 
low energies. As A grows a gap appears in 2D. Yet, if /i is 
set to emulate 3D, the DOS remains similar to the DA. As 
n c grows, the region over which the DOS spreads becomes 
a smaller fraction of A, since the number of zeroes grows 
with n c and the effect of the off-diagonal perturbations 
decreases, likely due to large cancellations in (El, In Fig. 
1 we show the DOS for n c = 40 and N = IOBq for the 
values of A = 0.01,0.5,4. The energy is rescaled by u c . 
The nearly diagonal case, A = 0.01, and A = 0.5 appear 



very similar. They both show a linear dependence in the 
DOS indicating qualitatively same behavior. The case 
A = 4 is qualitatively different. A gap is appearing at 
low energies and the next band has a small gap to the 
lowest one indicating the band-crossing is imminent. We 
can see that as A grows some structure is detected in the 
shape of the DOS that clearly develops for the A = 4 
case. These results demonstrate that within our pocket 
(H > H*,T < T*) the DOS is qualitatively similar to 
the DA and gapless and nearly gapless points {qf ,q"™} 
still dominate low energy properties. 

Our results have direct and observable consequences 
for a broad range of thermodynamic and dynamic prop- 
erties. Specific heat, ultrasound attenuation, STM and 
other forms of tunneling, and many other phenomena, 
will exhibit power-law-like behaviors, characteristic of 
low-energy BCS quasiparticles. All these diagnostic tools 
can now be used to explore the "extreme" pocket, identify 
gapless behavior and test our prediction that the high- 
and the low-field regime of type-II superconductors differ 
in this essential respect. 
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FIG. 1. DOS (in arbitrary units) for several values of A 
(see text). The inset shows the "extreme" pocket (shaded 
area) . 
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